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Resumen 


Las poblaciones de las especies invasoras pueden diferir en las características que determinan la 
intensidad de sus efectos sobre los ecosistemas nativos. El estudio de la variación morfológica aporta 
información valiosa sobre diferentes procesos evolutivos y ecológicos en las especies. Con el objetivo de 
evaluar los patrones de variación morfológica del caracol gigante africano (Achatina fulica), una especie 
invasora presente en el departamento del Valle del Cauca, se colectaron individuos en las zonas urbanas 
de 4 municipios (Buenaventura, Cali, Tuluá y Cartago). A estos especímenes se les realizó un análisis 
de morfometría geométrica y clasificación con base en el patrón de coloración. Se encontró que existe 
variación morfológica significativa dentro de la población, además del patrón de bandas o coloración, 
con un total de 13 patrones en toda la muestra. Se proponen las condiciones ambientales de cada 
localidad y los controles realizados por la autoridad ambiental como posibles causantes de la variación. 


Palabras clave. Estructura etaria. Invasión biológica. Moluscos. Morfología. Patrón de coloración. 


Abstract 


Populations of invasive species may differ in the characteristics that determine the intensity of their 
effects on native ecosystems. The study of morphological variation provides valuable information on 
different evolutionary and ecological processes. In order to evaluate the patterns of morphological 
variation of the African Snail (Achatina fulica), an invasive species present in the department of Valle 
del Cauca, individuals were collected in urban areas of 4 municipalities (Buenaventura, Cali, Tuluá and 
Cartago). An analysis of geometric morphometry and classification based on the pattern of bands was 
performed on the collected specimens. A significant morphological variation within the population 
was found, along with a total of 13 patterns of bands or coloration in the entire sample. Environmental 
conditions of each locality and controls exerted by the environmental authority are proposed as possible 
causes of variation. 
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Introducción 


Las especies invasoras son consideradas como una 
de las principales amenazas para la conservación 
de la biodiversidad global debido a que ejercen 
presiones sobre la fauna local a diferentes niveles, 
lo que puede desencadenar procesos de extinción 
local, regional o global (Williamson, 1996; Mack et 
al., 2000). El estudio de la variación morfológica de 
estas poblaciones aporta información valiosa sobre 
procesos evolutivos y ecológicos que se desarrollan 
a nivel intraespecífico e interespecífico. Incluso, ha 
sido propuesto que las variaciones morfológicas en 
estas poblaciones son el resultado de la forma en 
cómo se relacionan los individuos que la conforman 
con el ambiente, siendo modulada la respuesta 
por la información genética de los individuos 
(Bookstein, 1996). En este sentido, Parker et al. 
(2003) han propuesto que las variaciones espaciales 
y temporales en las características morfológicas 
de las especies invasoras son generadas por 
plasticidad o adaptación local, por lo que se hace 
necesario que estas variaciones sean consideradas 
para entender y manejar esta presión biológica 
externa en los sistemas afectados. 


El caracol gigante africano (Achatina fulica) es 
considerado una de las 100 plagas más peligrosas 
a nivel mundial por la Unión Internacional para 
la Conservación de la Naturaleza (UICN) debido 
a la alta proliferación, polifagia y simpatía por 
parte del ser humano (Patiño-Montoya y Giraldo, 
2017). La invasión de esta especie puede generar 
consecuencias devastadoras en el medio ambiente, 
como lo es el desplazamiento de especies nativas, 
ya que este caracol es un competidor fuerte 
y disminuye cultivos o plantas nativas por la 
polifagia. Adicionalmente, esta especie puede 
afectar la salud pública al ser intermediario de 
varios parásitos, como los nemátodos del genero 
Angiostrongylus, que producen enfermedades 
letales como la meningitis eosinofílica y la 
angiostrongilyasis abdominal (Alburquerque et al., 
2008, 2009; Liu et al., 2011). 


A. fulica es una especie originaria de Africa oriental, 
que ha venido expandiendo su distribución desde 
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hace aproximadamente 200 años (Fontanilla et 
al., 2014), y se ha establecido con éxito en todo el 
trópico y subtrópico del planeta (Alburquerque et 
al., 2008; Okon et al., 2012; De la Ossa-Lacayo et 
al., 2012; Avendaño y Linares, 2015). Como en la 
mayoría de moluscos, las diferencias en el tamaño, 
la forma y el patrón de bandas generalmente 
son atribuidas a respuestas adaptativas a nivel 
fisiológico y a la variación genética de la población, 
como una estrategia para afrontar las condiciones 
ambientales a las que se enfrentan (Okon et al., 
2012). Particularmente, para el caracol gigante 
africano ha sido reportada una alta variación en la 
forma de la concha en lugares donde la invasión 
es antigua (mayor a 20 años), ya sea de tipo 
primario o secundario (reintroducción o desde un 
sitio introducido) (Pattarmanon, 2004; Sobrepeña 
y Demayo, 2014a). Sin embargo, aún no se ha 
evaluado la variación en la forma de la concha 
de esta especie en lugares donde la invasión es 
reciente (menor a 20 años) como en Ecuador, Perú, 
Argentina, Venezuela o Colombia (Alburquerque 
et al., 2008; Alburquerque et al., 2009; Goldyn et 
al., 2016; Goldyn et al. 2017). En este trabajo de 
investigación se presentan los resultados del 
análisis cualitativo y cuantitativo de la variación 
en la forma y patrón de coloración de la concha 
del caracol gigante africano, en una región de 
reciente invasión, el departamento del Valle del 
Cauca, suroccidente de Colombia. 


Materiales y métodos 


Área de estudio. La recolección de individuos 
se realizó en las ciudades de Buenaventura 
(3%53"N-77" 050), Cali (3°26’N-76°31’O), Tuluá 
(4°05’N-76°12’0) y Cartago (4%44'N-75%54”) en 
el departamento del Valle del Cauca, Colombia. 
Buenaventura está situada en la región del Chocó 
Biogeográfico, a una altitud promedio sobre el 
nivel del mar de 7 m. La vegetación dominante 
es de tipo de bosque muy lluvioso tropical, con 
registro de precipitación anual promedio de 
7650 mm, humedad relativa promedio del 89 % 
y temperatura ambiental que varía entre 25 *C 
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y 28 °C (Rangel y Edgar, 2003). Las otras tres 
localidades se ubican en el valle geográfico del rio 
Cauca, a una altitud promedio de 1000 m s.n.m. 
La vegetación dominante es de tipo de bosque 
seco tropical, con registros de precipitación anual 
promedio de 900 mm, humedad relativa promedio 
del 70 % y temperatura ambiental que varía entre 
23 °C y 24 °C (Álvarez-López et al., 1984; Alberico 
et al., 2004) (Figura 1). 


Métodos de campo. Los individuos utilizados 
para el desarrollo de esta investigación fueron 
recolectados en el marco de un estudio de 
dinámica poblacional de caracol gigante africano 
en el Valle del Cauca (Giraldo et al., 2014). En 
términos generales, el área urbana de cada una de 
las localidades fue dividida en celdas de 2 km’. Al 
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Figura 1. Localización geográfica de las 


localidades muestreadas en el departamento 
del Valle del Cauca (Colombia): Buenaventura, 
Cartago, Tuluá, Cali. 
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interior de cada celda, se establecieron de manera 
aleatoria 3 parcelas de 16 m?, en donde se realizó 
el esfuerzo de búsqueda y recolección manual 
de caracoles vivos o conchas vacías, siguiendo 
las recomendaciones de seguridad para la 
manipulación de esta especie definidas por la 
autoridad ambiental. 


Análisis morfométrico y de patrón de bandas. 
Con el propósito de establecer la estructura de 
talla de la población de A. fulica en la zona de 
estudio, se midió la altura total de la concha de 
cada individuo capturado y se construyó un 
histograma de talla considerando intervalos de 
10 mm. Se desarrolló un análisis de morfometría 
geométrica para evaluar la variación de la forma 
utilizando individuos adultos jóvenes con tallas 
entre 40 a 60 mm (Simáo y Fischer, 2004), con el 
propósito de reducir la pérdida de información por 
desgaste de la concha y el efecto del crecimiento 
alométrico que presentan este tipo de especies. 


Se realizó un registro fotográfico digital desde la 
vista ventral de 47 individuos: 12 de Buenaventura, 
11 de Cali, 12 de Tuluá y 11 de Cartago y se 
utilizaron las herramientas de procesamiento 
digital para análisis morfométrico Thin-Plate 
Spline (TPS)  (http://life.bio.sunysb.edu/ee/ 
rohlf /software.html). Se definieron 10 puntos 
de referencia distribuidos en el ápice, la sutura, 
la columela y la apertura de la concha (Figura 2). 
Para la transformación de las imágenes a formato 
TPS se utilizó TPSultil®. Posteriormente, se 
estableció la posición de los puntos de referencia 
para cada individuo utilizando TPSDig2® y 
se procedió a realizar un análisis general de las 
deformaciones relativas utilizando la rutina 
TPSrelw®. Este análisis se fundamenta en la 
comparación de la posición de los puntos de 
referencia, de todos los individuos, respecto a 
una forma consenso (Toro et al., 2010). Para 
determinar si la forma de las conchas es 
independiente del tamaño, se realizó un análisis 
de regresión con los dos primeros componentes 
resultantes del Análisis de Componentes Princi- 
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Figura 2. Landmarks digitali- 
zados en orden numérico sobre 
las estructuras representativas 
de la concha de A. fulica: ápice, 
espiras, vuelta del cuerpo y 
apertura. 


pales (ACP) de las deformaciones relativas como 
variables dependientes y el tamaño del centroide 
como variable independiente. 


Para evaluar la variación entre localidades se 
realizó un análisis discriminante con base en 
las deformaciones relativas y posteriormente se 
compararon las localidades a través de un análisis 
de varianza (ANOVA) utilizando los valores 
de las dos primeras funciones discriminantes. 
Finalmente, se utilizó TPSsplin® para visualizar 
los cambios relativos de la forma de la concha que 
más contribuyeron a las funciones discriminantes 
establecidas en el análisis previo. Los análisis 
estadísticos se realizaron con los paquetes Vegan 
y MASS (Venables y Ripley, 2002). 


Para el análisis del patrón de bandas se siguió 
la clasificación de Sobrepeña y Demayo (2014a), 
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quienes definieron 14 patrones de bandas para 
la especie. Cada individuo fue clasificado de 
acuerdo con el patrón de bandas, describiéndose 
el patrón predominante para cada localidad 
y en toda la zona de estudio. Con el fin de 
evaluar si existe relación entre los patrones de 
bandas y las localidades, se realizó un análisis 
de contingencia y posteriormente pruebas chi- 
cuadrado secuenciadas en el programa PASTO, 
para establecer cuales localidades eran la fuente 
de variación. 


Resultados 


Análisis morfométrico. Para realizar el análisis 
de variación de talla se procesaron un total de 192 
conchas de A. fulica: 24 de Buenaventura, 39 de 
Cali, 80 de Cartago y 49 de Tuluá. Las poblaciones 
de Buenaventura y Tuluá presentaron los mayores 
promedios de altura total con 63,85 + 8,83 mm y 
63,25 + 12,97 mm respectivamente, seguidos de 
Cali con 53,71 + 12,97 mm y Cartago con 39,18 
+ 11,02 mm. De acuerdo con los registros de 
tamaño, los adultos jóvenes de la zona de estudio 
no superaron el 60 % del total de individuos 
capturados (Figura 3). 


Los dos primeros componentes del ACP 
explicaron el 65 % de la varianza contenida en las 
deformaciones relativas. La regresión entre estos 
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Figura 3. Estructura etaria según Simao y Fisher 
(2004) de las conchas de Achatina fulica evaluadas 
en cada localidad. 
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componentes y el tamaño del centroide no fue 
significativa (F,,) = 0,37,p = 0,68), lo que indica 
que no existió una variación significativa en la 
forma de la concha en el rango de tamaño de los 
individuos estudiados (40-60 mm), por lo que las 
variaciones en forma que sean detectadas serían 
independientes del tamaño. 


El análisis discriminante indicó que el 70 % de la 
varianza en la forma de la concha de A. fulica estuvo 
contenida en las dos primeras dimensiones. Las 
deformaciones relativas (RW) que más aportaron 
a la primera función discriminante (FD1) fueron 
el 2X y el 2Y, las cuales representaron el ancho 
de la apertura amplia, y el ancho y el tamaño de 
la espira. La segunda función (FD2) representó 
el tamaño de la espira y el ancho de la vuelta 
del cuerpo angosta (RW 4Y y 6Y). Por último, la 
tercera función (FD3) la orientación de la espira y 
del ápice (RW 1Y-2Y) (Figura 4, Tabla 1). 


La primera función discriminante (FD1) 
proporcionó evidencia significativa de la 
diferencia de la forma de la concha entre las 
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< 0,05), permitiendo identificar una diferencia de 
la forma entre los especimenes de Buenaventura 
con los especimenes de las localidades del valle 
geografico del rio Cauca (Figura 3). Entre estas 
últimas localidades la FD1 no sugirió la presencia 
de diferencias en la forma de la concha (p = 0,88 
y p = 0,08). Sin embargo, al incluir en el análisis 
la variación de la segunda (FD2, F.,,, = 8,08, 
p 0,0001) y tercera (FDS, F „44 - 8,08, p = 0,0001) 
función discriminante, se detectaron diferencias 
significativas en la forma de la concha de A. fulica 
entre Tuluá y Cartago (p < 0,05) y Tuluá con Cali 
(p < 0,05). 


En términos generales, la forma de la concha de los 
individuos de A. fulica que habitan en la ciudad de 
Buenaventura tiende a tener una apertura amplia 
y una espira delgada y alargada, lo que se traduce 
en una concha más estilizada con respecto a las 
poblaciones del valle geográfico del río Cauca. Por 
otra parte, la forma de la concha de las poblaciones 
de A. fulica que habitan en el valle geográfico del 
río Cauca tienden a variar entre sí en el ancho de 
la vuelta del cuerpo y en la orientación del ápice 





diferentes localidades de estudio (F,,, - 17,14, p (Figura 5). 
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Figura 4. Análisis discriminante para las poblaciones de A. fulica muestreadas en 


cuatro municipios del Valle del Cauca. 
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Tabla 1. Aportes de las deformaciones relativas 
a las funciones discriminantes FD1, FD2 Y 
FD3. LD1, LD2 Y LD3: Coefficients of Lineal 
Discriminants. 















































LD1 LD2 LD3 
1X -8,35E-13 6,97E-13 -1,14E-12 
1Y -4,20E-13 -3,94E-13 2,66E-12 
2X 1,02E-12 -5,92E-13 -1,04E-12 
2Y -1,09E-12 3,26E-13 1,54E-12 
3X -6,51E-13 1,13E-13 -3,68E-13 
3Y 1,80E-13 8,14E-13 5,44E-13 
4X 2,54E-13 -1,53E-13 3,59E-13 
4Y -6,89E-13 -2,68E-12 1,33E-13 
5X 2,71E-13 -5,68E-13 -1,19E-12 
5Y -4,74E-13 2,01E-13 -5,39E-13 
6X 7,55E-13 1,40E-13 2,78E-13 
6Y 4,83E-14 1,75E-12 -3,54E-13 
7X -1,04E-12 -5,31E-13 7,94E-13 
7Y -1,01E-13 -1,96E-13 5,23E-14 








LA AA 
Figura 5. Formas de la concha de A. fulica 
asociadas a las funciones discriminantes 
FD1, FD2, FD3. A) PW1, B) PW2, C) PW4, 
D) PW6. 
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Patrón de coloración. De los 14 patrones de 
coloración que han sido descritos para el caracol 
gigante africano por Sobrepeña y Demayo (2014a), 
el único que no se encontró en la zona de estudio 
fue el patrón # 14. Las localidades de Cali y 
Buenaventura obtuvieron la menor cantidad de 
patrones diferentes (6), a diferencia de Cartago y 
Tuluá en donde se registraron 11 y 10 patrones de 
bandas, respectivamente. El patrón #5 fue el más 
abundante (25,5 %) en el área de estudio, seguido 
del patrón #1 (22,9 %). Sin embargo, este último 
patrón fue el más abundante para las localidades 
de Cartago y Tuluá (> 30 %), mientras que el patrón 
#5 fue el más abundante en Cali (33,3 %) y #7 (29,2 
%) y #9 (20,8 %) fueron los más abundantes en 
Buenaventura (Tabla 2, Figura 6). Es importante 
destacar que en ninguna de las localidades de 
estudio se registró el total de los patrones de 
bandas registradas en todo el Valle del Cauca (13 
patrones). 





Figura 6. Patrones de coloración según la 
clasificación de Sobrepeña y Demayo (2014a) 
predominantes en los cuatro municipios del 
Valle del Cauca. A) Patrón 5 (Cali), B) Patrón 7 
(Buenaventura) y C) Patrón 1 (Cartago y Tuluá). 
No se estableció independencia entre el patrón de 
coloración y las localidades (x?¿¿, = 75,7, p 0,0001), 
siendo la frecuencias de los patrones de coloración 
significativamente diferentes entre localidades 
Wes = 30, p valor >> 0,05). Este resultado sugiere 
que cada localidad tuvo un patrón específico 
de banda dominante y además cada localidad 
exhibió una frecuencia diferente de patrones de 
bandas presentes. 
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Tabla 2. Frecuencia expresada en porcentaje de los patrones de bandas 
presentes en las conchas de Achatina fulica de las localidades de Buenaventura, 
Cali, Tuluá, Cartago y el total. n: número de conchas. 















































Buenaventura Cali Cartago Tuluá Total 

Patrón 1 0,0 2,6 ss 36,7 AS 
Patrón 2 4,2 0,0 3,8 0,0 2,1 
Patrón 3 16,7 10,3 3,8 8,2 7,8 
Patron 4 0,0 0,0 2,5 6,1 2,6 
Patrón 5 125 399 25,0 26,5 255 
Patrón 6 0,0 0,0 0,0 2,0 0,5 
Patron 7 292 15,4 6,3 0,0 9,4 
Patrón 8 16,7 12,8 11,3 2,0 9,9 
Patron 9 20,8 25,6 15 10,2 118,5 
Patrón 10 0,0 0,0 2,5 4,1 211 
Patrón 11 0,0 0,0 0,0 220 0,5 
Patron 12 0,0 0,0 5,0 2,0 2,6 
Patron 13 0,0 0,0 is 0,0 0,5 

n 24 39 80 49 192 





Discusion 

Entender cómo reaccionan las especies invasoras, 
como el caracol gigante africano, a los factores 
bióticos y abióticos en un entorno foráneo es 
importante para establecer planes de control 
y manejo de este tipo de especies (Sakai et al., 
2001; Davidson et al., 2011; Avendaño y Linares, 
2015). Particularmente para los moluscos, las 
variables que más aportan información sobre las 
interacciones con el ambiente son las variaciones 
en tamaño y forma, de ahí la importancia de 
conocer los patrones de variación de estos atributos 
morfológicos en un área de interés (Mack et al., 
2000; Davidson et al., 2011). 


En este trabajo no se encontró evidencia en 
la fracción de tamaño estudiada de tendencia 
alométricas, por lo que es de esperar independen- 
cia entre la forma y el tamaño. Bajo este escenario, 
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los patrones de variación en la forma y en colo- 
ración estarían siendo modulados por factores 
diferentes al tamaño. Factores como las condicio- 
nes climáticas, la estructura y distribución 
espacial de zonas verdes o lotes baldíos de 
las ciudades, múltiples eventos de introducción 
e incluso las medidas de control que realizan 
las autoridades ambiental locales y regionales, 
estarían afectando de manera significativa el 
tamaño y la forma de estos moluscos (Raut y 
Baker, 2002; Alburquerque et al., 2009; Avendaño 
y Linares, 2015). Además, las variaciones espacia- 
les en la longitud de la concha encontradas en 
este estudio pueden deberse a las diferencias 
climá-ticas entre los ambientes ocupados por el 
caracol gigante africano en el Valle del Cauca, 
que va de un ambiente húmedo del Chocó 
Biogeográfico (Buenaventura) a un ambiente 
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seco del valle geográfico del río Cauca (Cali, 
Tuluá, Cartago). Este resultado es similar al 
reportado por Giraldo et al. (2014) y Avendaño 
y Linares (2015), quienes proponen que las 
características ambientales presentes en la ciudad 
de Buenaventura tenderían a favorecer la tasa 
de crecimiento de esta especie, ya que no está 
sujeto a un fuerte estrés hídrico, lo que resultaría 
en conchas más grandes, con apertura amplia y 
espira delgada y alargada. 


La disponibilidad de microhábitats que puede 
ocupar el caracol gigante africano en un ambiente 
urbano (lotes baldíos, jardines, fuentes hídricas 
o parques) puede ser un elemento que estaría 
promoviendo la variación del patrón de bandas al 
interior de cada localidad. De acuerdo con Giral- 
do et al. (2014), la disponibilidad de microhá- 
bitats favorables para el caracol gigante africano 
en las localidades de estudio es significativamente 
diferente, siendo mayor la disponibilidad de 
lotes baldíos en Buenaventura, mientras que en 
Cali, Tuluá y Cartago los jardines, separadores 
viales y zonas verdes asociadas a parques fueron 
más frecuentes. Como se ha reportado en otros 
moluscos, la complejidad estructural del ambiente 
donde se desenvuelvan está relacionado con la 
forma de crecimiento (Flores-Campaña et al., 
2012). Por lo tanto, en el caso del caracol gigante 
africano, su relación con la disponibilidad de 
microhábitats en los centros urbanos con procesos 
de invasión reciente podría estar modulando el 
tamaño y el ancho de la espira, así como la vuelta 
del cuerpo (Tomiyama, 1992; Pattarmanon, 2004; 
Alburquerque et al., 2009; Sobrepeña y Demayo, 
2014b). 


Debido a la naturaleza invasora del caracol 
gigante africano, es probable que al departamento 
del Valle del Cauca estén ingresando individuos 
de diferentes lugares, en diferentes temporadas 
e incluso con diferente intensidad (Fontanilla et 
al., 2014; Patiño-Montoya y Giraldo, 2017). Lo 
anterior, sumado al control diferencial sobre 
clases específicas de tamaño de la población 
(adultos) como consecuencia de las acciones 
de control implementada por la autoridad 
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ambiental, explicaría la constante movilización 
de los individuos dentro de la población que se 
encontró en el análisis discriminante. Finalmente, 
el patrón de coloración descrito en este trabajo 
para el caracol gigante africano en el Valle del 
Cauca podría ser el resultado de la interacción 
conjunta de forzantes genéticos y ambientales que 
estarían modulando las frecuencias de ocurrencia 
de patrones de bandas en cada localidad, e 
igualmente su variabilidad espacial (Sobrepeña 
y Demayo, 2014a, 2014b; Mack et al., 2000; 
Alburquerque et al., 2008; Pattarmanon, 2004; De 
la Rosa et al., 2010). 


Conclusión 


El tamaño del ápice, de la apertura y la forma 
de la vuelta de la concha del caracol gigante 
africano (Achatina fulica) en el Valle del Cauca 
fueron las fuentes de variación de la forma 
entre las localidades de estudio, además de que 
se establecieron diferencias espaciales signifi- 
cativas en los patrones de coloración de la con- 
cha. Considerando los resultados obtenidos, 
probablemente el efecto de esta especie invasora 
sobre los hábitats urbanos del Valle del Cauca 
no es uniforme. Por lo tanto, la capacidad de 
resiliencia de los sistemas urbanos afectados por 
la invasión del caracol gigante africano también 
sería variable. 
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